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SUMMARY 

Global solar radiation measurements were available as chart 
recordings at two stations, Frelighsburg and L'Acadie. These recordings 
were converted to daily radiation values using a digitalizing table and 
conversion constants appropriate to each station. Three methods of 
estimating missing data were tested. The first method consisted of 
correlating global solar radiation from a neighbouring station, 
Jean-Brébeuf, with measurements available at Frelighsburg and L'Acadie. 
At Jean-Bréboeuf, solar radiation was available since 1964 on a daily 
basis from Environment Canada, and correlation with this station 
provided the best estimates of missing data at Frelighsburg and 
L'Acadie. However, if reliable data are not available at a neighbouring 
station, the third method tested, based on relating hours of sunshine to 
global radiation, also provided good estimates of missing data. The 
second method involved relating air temperature to global radiation, but 


this method is not recommended for use during the winter months. 
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RESUME 

Les mesures de rayonnement global sont disponibles sous forme de 
tracés sur carte à deux stations, Frelighsburg et L'Acadie. Ces tracés 
ont été convertis en valeur énergétique journalière à l'aide d'une table 
de digitalisation en utilisant une constante de calibration propre à 
chaque appareil. Trois méthodes ont été évaluées afin d'estimater les 
données manquantes. La première méthode consiste à établir des 
corrélations entre des stations voisines. En plus de Frelighsburg et 
L'Acadie nous avons utilisé les données de rayonnement global à 
Jean-Brébeuf. La corrélation provenant des données de cette station, 
disponibles depuis 1964 à Environnement Canada, a fourni 1a meilleure 
estimation des données manquantes à Frelighsburg et L'Acadie. 
Cependant, si des données de qualité ne sont pas disponibles à une 
station voisine, la troisième méthode évaluée, basée sur la corrélation 
entre l'ensoleillement et le rayonnement global, fournit aussi une bonne 
estimation des données manquantes. La seconde méthode consistant à 


utiliser les données de température de l'air et de rayonnement global, 
ÿ 8 


n'est pas recommandée durant les mois d'hiver. 
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Introduction 

Le rayonnement solaire global est tout le rayonnement solaire 
descendant de courtes longueurs d'ondes, direct et diffus, reçu sur 
une surface horizontale plate (Service de l'environnement 
atmosphérique, 19/78). En recherche agrométéorologique, ces mesures 
ne sont pas abondantes. Et pourtant, elles ajoutent une information 
indispensable sur l'énergie disponible à la photosynthèse et la 
transpiration. Aussi était-il intéressant de pouvoir ajouter à 
notre banque actuelle les données provenant de L'Acadie et de 
Frelighsburg. 

Ces données étaient cependant encore sous leur forme originale 
c'est-à-dire des courbes sur des cartes. Elles ont donc été 
intégrées à l'aide d'une table de digitalisation puis converties en 
TR En effet, ces unités sont celles employées par le Service de 
l'environnement atmosphérique depuis 1978 (Service de 
l'environnement atmosphérique, 19/78). Bien que dans les articles 
cités dans ce bulletin, les unités employées soient généralement les 
nee (Rare = 0,04186 SH LA à les Monte sont les unités 
conformes au système international (SI). 

Ce travail a été ensuite complété par une évaluation de trois 
méthodes mathématiques pour l'estimation du rayonnement global lors 


des journées manquantes. 


Evaluation de la constante de calibration 


2.1 Méthodologie 
À Frelighsburg, l'enregistreur fonctionne depuis 1979. Il 


s'agit d'un pyranographe mécanique modèle 4001 fabriqué par 


Sierra-Misco, Inc (Californie, E.-U.) . L'instrument utilisé à 
L'Acadie est un pyrhéliographe fabriqué par Belfort (Baltimore, 
E.-U.).. Il a fonctionné de 1973 à 1980. Depuis, il-est- hors 
d'usage. 

La méthode pour convertir un tracé sur carte en valeur 
énergétique journalière est simple. On mesure d'abord la surface 
sous la courbe du tracé sur carte à l'aide d'une table de 
digitalisation. Cette intégration a été faite à l'aide de 
l'équipement informatique du Système d'informatique des sols au 
Canada (SISCan). 

Par la suite, ces valeurs de surface en a sont converties en 
Mn d en utilisant une constante de calibration propre à chaque 


instrument (section 2.2). 


2.2 Calcul de la constante de calibration 
2.2.1 Frelighsburg 
Le manuel d'utilisation du pyranographe mécanique 


(Sierra-Misco, 1978) nous fournit l'information suivante: 
Rayonnement solaire = K * Lecture sur la carte 


ae 0,539Mg=cal 


2 
cm -min 
Le terme “lecture sur la carte” refère à la hauteur de 
la courbe à un moment donné. Ceci signifie que pour une 
hauteur de 1 cm sur la carte de rayonnement, le sol reçoit 


0,539 eallen tn 


=3%— 


Lors de l'intégration des courbes avec la table de 
digitalisation, nous obtenons une mesure de surface par jour. 


Or, 
Surface = Hauteur * longueur 


où longueur refère au déplacement du baril soit 4 cm par jour. 
Donc pour obtenir la valeur de rayonnement global pour 


une journée, nous calculons: 
À 2 
Rayonnement solaire (MJm }) = K * Hauteur 


Surface 


RE ———_—]—— 
longueur 


Kg £-cal Ce 00,05196 Mon x Surface _. 


cm -min cm 1 nl en & cm 
—: 
CN RE Emi, 1 | 0,04186 Mim 
cm -min DE hr cm lg-cal-em 


2 
surface cm 
& cm 


2 360 


0,02256 MJm “ * _ x Surface 


Rayonnement solaire = (8,123 * surface) Mn 


où surface est la valeur fournie par 1a table de 


digitalisation. 


2.2.2 L'Acadie 

Le fabricant du pyrhéliographe de L'Acadie (Belfort, 
1965) a fourni une méthode différente de celle utilisée à 
Frelighsburg pour le calcul de la constante. 

On calcule d'abord la valeur énergétique de la surface 
maximale pour une journée. Dans ce cas, / cm de hauteur sur la 
carte équivaut à 3 ere ue Après conversion en valeur 


journalière, nous obtenons que / cm de hauteur équivaut à: 


2 


3 g-cal _ 3g-cal , 1440 min 0,04186 MJm_ 


* 
2 2-1 ; 7 
cm -min cm -min JE g-cal-cm 


2 


180,8352 MJm. 


Or si on intègre la surface maximale à l'aide de la table de 


digitalisation, on obtient: 


Surface = Hauteur * longueur 


% 4,5 cm 


Il 
= 
Q 
5 
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Finalement, la constante K, se déduit ainsi: 


K, =1180;;:035 unie Rene 5,7408 MIE par _. sur la carte. 


À 31P5 ion 


Pour connaître la valeur de rayonnement global sous la surface 


réelle 


st 


K, * surface 


Rayonnement global à 


Donc, 


5,7408 * surface AE 


Rayonnement global 


où surface est la valeur obtenue par la table de 


digitalisation. 


3. Estimation des données manquantes 
3.1 Méthodologie 

Il arrive fréquemment qu'avec le type d'instrument employé, il 
manque une ou plusieurs journées de données. Trois méthodes ont 
été évaluées pour l'estimation des données manquantes. 

La première méthode consiste à établir des corrélations entre 
des stations voisines. En plus de Frelighsburg et L'Acadie, nous 
avons utilisé les données de rayonnement global à Jean-Brébeuf. Ces 
données sont disponibles depuis 1964 et ont été obtenues 
d'Environnement Canada. Cette station a été retenue pour la qualité 
de ses données et pour sa proximité des sites à l'étude. 


Les équations sont de la forme: 


YE=N a FI Dx 


où a est l'ordonnée à l'origine en Mme 
b est la pente 
—2 
x le rayonnement global à la station observée (MJm ) 


Y le rayonnement global à la station estimée (nee) 


La deuxième méthode vérifie l'approche exponentielle suggérée 


par Richardson (1984). Cette approche s'écrit: 


Q 


_S = À x (TMX-TMN)P 
QE 


où TMX = température maximale quotidienne en degrés Celsius 


TMN = température minimale quotidienne en degré Celsius 


2 


QE = rayonnement global quotidien (MJm ‘) 
Q, = rayonnement globale théorique au sommet de l'atmosphère 
(He) 


A,B = coefficients 


Des variantes à cette équation ont aussi été évaluées telles que 


A * (TMX-TMN)? 


Q£ = 
0570 = À + B * (TMX-TMN) 
QE AUS (TMX-TMN ) 


Les valeurs de température et rayonnement utilisées sont celles 
de Frelighsburg et Jean-Brébeuf. 

La troisième méthode, dont une revue de littérature est 
présentée par Hayhoe (1980), utilise la relation linéaire entre le 


rayonnement global et l'ensoleillement: 


a/Q = À + B (n/N) 


où n = nombre d'heures réelles d'ensoleillement 
N = nombre d'heures théoriques d'ensoleillement 
À = ordonnée à l'origine 


B = pente de la droite 


Les seules données utilisées sont celles de Jean-Brébeuf. 


3.2 Résultats et discussion 
3.2.1 Corrélations simples entre deux stations 
Les résultats des diverses corrélations apparaîssent au 
tableau 1. Deux séries de corrélation ont été faites en 
prenant les stations par paires. La première série utilise la 
même période (octobre 1979 à décembre 1980) pour toutes les 
paires; la seconde série utilise toutes les données disponibles 
pour une paire de station. Les figures 1 à 5 illustrent la 
relation qui existe entre chaque paire. Sur ces figures, “cm ” 
réfère à la surface sous la courbe telle que mesurée avec la 
table de digitalisation. 

Si on considère la même période d'observation (281 
données) et en se basant sur les corrélations, il est 
préférable d'utiliser les données de Jean-Brébeuf si l'on veut 
estimer le rayonnement à la station de L'Acadie ou à celle de 
Frelighsburg (différences significatives àa =0,05). Cependant, 
la corrélation entre deux stations diminue si on augmente la 
période couverte de sorte que les différences entre les 


corrélations ne sont plus significatives (seuil a = 0,01). Les 
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données de L'Acadie (ou Frelighsburg) sont alors aussi valables 
que celles de Jean-Brébeuf pour estimer Frelighsburg (ou 


L'Acadie). 


3.2.2 Approche exponentielle et ses variantes 

Le tableau 2a montre les résultats obtenus pour les 
diverses approches testées à Frelighsburg. L'utilisation de la 
fraction Q/Q plutôt que de la radiation seule dans les 
équations exponentielles et linéaires, n'améliore pas de façon 
significative (seuil a=0,05) l'estimation des valeurs de 
rayonnement global à la station. La forme exponentielle donne 
des corrélations significativement plus élevées (a=0,05) que la 
forme linéaire sauf aux mois d'octobre et novembre. 

L'effet des mois est nettement le facteur prépondérant, 
les mois d'été présentant des corrélations plus élevés que les 
mois d'hiver. Notons que la corrélation trouvée en janvier 
dans les variantes linéaires, n'est pas significative (seuil 
a=0 ,05),. 

Les calculs faits avec les données de Jean-Brébeuf 
(tableau 2b) conduisent aux mêmes conclusions. Dans ce cas 
cependant, il n'y a pas de différence significative (seuil 
a=0,05) entre l'approche linéaire et exponentielle. L'effet 
des mois ressort ici aussi: en effet, les mois de janvier, 
février et décembre ne montrent pas de corrélation 


significative (seuil a=0,05) entre le rayonnement estimé par 


12 
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SAS 


les températures et le rayonnement global observé et ce, dans 


les quatre équations testées. 


3.2.3 Relation avec l'ensoleillement 

Les résultats obtenus sont présentés au tableau 3. Les 
corrélations s'avérent toutes élevées et significatives au 
seuil a=0,05. Elles sont du même ordre de grandeur que celles 
calculées par la première méthode au tableau 1. La pente et 
l'ordonnée à l'origine pour l'ensemble de l'année se comparent 
aux valeurs trouvées par Hayhoe (1980) pour quatre stations 


canadiennes. 


Conclusion 


Les tracés de rayonnement global à Frelighsburg et L'Acadie, 
ont été intégrés et sont maintenant disponibles en valeurs 
énergétiques (Mini) sur ruban magnétique à la section 
d'agrométéorologie, Agriculture-Canada. 

Le nombre de données manquantes nous a incité à tester diverses 
méthodes pour les estimer. La meilleure approche consiste à estimer 
la valeur manquante à partir de la donnée d'une station voisine. 
Cependant, si aucune donnée n'est ainsi disponible, l'utilisation de 
l'ensoleillement est la seconde solution. Finalement, une troisième 
solution consistant à utiliser les données de température maximale 
et minimale quotidienne peut être employée en dernier recours. 
Cependant, cette dernière approche n'est pas recommandée durant les 


mois de décembre, janvier et février. 


s'TRE 


Tableau 3. Relation entre la radiation globale et l'ensoleillement à 


Jean-Brébeuf pour la période 1964 à 1983 


Période 


Printemps 
Eté 


Automne 


Hiver 


Année 


Hayhoe (1980) 


Ordonnée à l'origine 


02223 
0,227 
05 232 


0,321 


0,258 


0523 


Pente 


0,538 
0,487 
0,632 


0,583 


0,525 


0,54 


Corrélation (%) 


93,7 
91,3 
90,3 


90,0 


909 


90,0 
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Annexe l 


Estimé des valeurs manquantes 


(Mme) 


Légende: 

(1) valeur observée à la station 

(2) valeur estimée à partir de la station voisine (méthode 1) 
5) valeur estimée à partir de Jean-Brébeuf (méthode 1) 


(3)*: valeur estimée à partir d'une autre méthode (méthode 3 ou 2) 


les 


Pour estimer les données manquantes, nous avons finalement retenu 
équations suivantes: 


1,/276 + 0,8148 * Jean-Brébeuf 
1,8901 + 0,8861 * L'Acadie 


il 


Frelighsburg 
Frelighsburg 


L'Acadie = 1.3468 + 0,8075 * Jean-Brébeuf 
L'Acadie 0,5385 + 0,8877 * Frelighsburg 


Lorsque les données de rayonnement sont manquantes aux trois stations, 
nous avons utilisé la méthode basée sur l'ensoleillement ou si cette 
donnée est aussi manquante, les températures quotidiennes à 
Frelighsburg. 
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Annexe 2 


Les données de rayonnement global sont disponibles sur ruban à la section 


d'agrométéorologie, Ottawa, Ruban DL4876, label 13, 
DSN=AG775. PMET. ANNEX3. RAPPORT. ACAFREL. BREBOEUF. 


Chaque enregistrement a une longueur de 7/3 caractères. 


Le format des données est le suivant: 


STATION DATE VALEUR SURFACE LONGUEUR FRELI/ J. -BREBEUF 


L'ACADIE 
A13 4212 -F10.3 F6.3 F6.3 F10:3 P10.3 
STATION: nom de la station: Frelighsburg ou L'Acadie 
DATE : date selon le format 19AA-MM-JJ 
VALEUR : valeur de rayonnement, MJm observé à la statio 
SURFACE : valeur obtenue de la table de digitalisation, cm 


LONGUEUR: 


AUTRE 


F10.3 


, cm (base) 


FRELI/L'ACADIE: valeur de rayonnement estimée à partir de l'autre station 


JEAN-—BREBEUF : valeur de rayonnement estimée à partir de Jean-Brébeuf 


AUTRES : valeur de rayonnement estimée par une autre des méthodes 


décrites dans le texte. 
I valeur entière 
A: caractère alphanumérique 
F valeur réelle avec décimale 


Lorsqu'une donnée est manquante, elle est identifiée par -9.999 


Annexe 3 


Programme pour le calcul des valeurs théoriques de Q0 et N 


Sous-routine qui calcule les valeurs théoriques de rayonnement global au 
sommet de l'atmosphère et de nombre d'heures maximales d'ensoleillement 
(Adaptée de Robertson et Russelo, 1967). 


SUBROUTINE RAYON (DAYLEN, Q0 ,ALAT, JD,JF) 


C—— DAYLEN = Nombre d'heures d'ensoleillement théorique 

ee Q0 = Rayonnement global au sommet de l'atmosphère 
(Ce ALAT = Latitude en degrés, centième 

(CS JD = Début de la période; 1=l janvier 

C— JF = Fin de la période de calcul; 365 = 31 décembre 


DIMENSION DAYLEN(1),Q0(1) 


PI = 3,141592654 
PITHE = 2.*P1/365. 
BLAT = ALAT *PI/180.0 
DAYF = 24.0/(2.0*PI) 
SOLC = 1.94%60. 


DO 2000 J=JD,JF 
THET=J*PITHE 
DEC=0. 3964E-00+0. 3631E+01*SIN(THET)-0. 2297E+02*COS(THET) 
+0.3838E-01*SIN(2.0*THET)-0. 3885E-00*CO0S( 2.0*THET) 
+0.7659E-01*SIN(3.0*THET)-0. 1587E-00*C0S(3.0*THET) 
-0,1021E-01*C0S(4.0*THET) 
DEC=DEC*PI /180. 
H=ACOS((-.01454 -SIN(BLAT)*SIN(DEC))/(COS(BLAT)*COS (DEC) ))*DAYF 
H1=ACOS ( (-S IN(BLAT ) *S IN(DEC) )/(COS(BLAT)*COS(DEC) )) 
RADV=1 .0-0.0009464*SIN(THET)-0.01671*COS(THET)-0. 0001489% 
COS(2.0*THET)-0.00002917%SIN(3.0*THET)-0. 00003438 
kCOS(4.O*THET) 
QO(J) = 2.0*SOLC *DAYF/(RADV*RADV)*(SIN(BLAT)*SIN(DEC)*H1 
+COS(BLAT)*COS(DEC)*SIN(H1))*0.04186 
DAYLEN(J) =2.0%H 
2000 CONTINUE 


RETURN 
END 


